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Abstract 

We study a family of classical solutions of modified sinh-Gordon equation, dzdzi] — e^^ + 
p{z)p{z) e~'^^ = with p{z) = z^" — s^". We show that certain connection coefficients 
for solutions of the associated linear problem coincide with the Q-function of the quantum 
sine-Gordon (a > 0) or sinh-Gordon (a < —1) models. 
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1 Introduction 

In over three decades of study of quantum integrable systems, a remarkable (and largely mysterious) 
relation to classical integrable equations was observed in a number of different contexts. The first 
such relation was discovered long ago by Barouch, Tracy, and McCoy pj, who have derived the spin- 
spin correlation function in the Ising Field Theory in terms of special Painleve III transcendent, 
i.e. the solution of the differential equation, 

g + ;^-i-i>(2/) = o (1.1) 

which decays at r — )• oo and is singular as — log(— rlogr) at r — )• 0. The derivation relies on the 
free-field nature of the Ising Field Theory, and it is still unknown if this powerful result can be 
generalized in any simple way to interacting integrable quantum field theories. 

However, relation to classical integrable equation has surfaced later in analysis of the dilute 
self-avoiding polymer problem [2l|3], which is related to quantum sine-Gordon model at special 
value of the coupling parameter. It was found that the off-critical partition function of a self- 
avoiding polymer loop on an infinite cylinder is expressed exactly through another solution of the 
same Painleve III equation p. ID . this time with the singularity — |logt + 0{1) at small t. This 
elegant result can be attributed to supersymmetry of the problem - it can be reformulated in terms 
of quantum sine-Gordon model at special value of its coupling constant = | in ()1.2p below), 
where it exhibits M = 2 supersymmetry. Indeed, the derivation in [4J has its roots in deep analysis 
of 2D N = 2 super symmetric field theories [5J. At the same time, the finite-size partition function 
is generally much simpler object as compared to the correlation functions. In particular, TBA 
technique [7] can be employed to determine this quantity in any integrable quantum field theory 
(supersymmetric or not) as long as its S-matrix is known [8]. Therefore one may expect that some 
more or less direct extension of the results of ^i3j to generic integrable theories is possible. 

In this paper we propose such extension to the case of the quantum sine-Gordon model 
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{dt^f-{d^vf + 2// cos (/3(^) (1.2) 
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at generic values of the coupling parameter < 1. Here sets the mass scale, ^ ~ [mass]^^^'^ 
We will consider the theory in finite-size geometry, with the spatial coordinate x va. ip = ip{x, t) 
compactified on a circle of a circumference R, with the periodic boundary conditions 

ip{x + R,t) = ip{x,t) . (1.3) 

Due to the periodicity of the potential term 2/i cos (/?</?) in ()1.2p in ip, the space of states Ti splits 
into orthogonal subspaces T-L^: characterized by the "quasi-momentum" k, 

(^^^ + 27r//3: | ) ^ e^^'M ^-t ) (1-4) 

for \ G Tij^. We call A;- vacuum the ground-state of the finite-size system (jl.2p in the sector T-Lk. 

The quantum field theory (jl.2p is integrable, in particular it has infinite set of commuting local 
Integrals of Motion (IM) l2n-i, ^2n-i, 2n = 2, 4, 6, . . . being the Lorentz spins of the associated 
local densities [6]. Of primary interest are the /c-vacuum eigenvalues l2n-iik \ R), l2n-i{k \ R), 
especially the fe- vacuum energy Ii + Ii. In principle, these quantities are accessible through the 
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Thermodynamic Bethe Ansatz (TBA) technique, but the most efficient approach is the Destri-De 
Vega (DDV) equation [9l[10j (Similar equation was earlier derived in the lattice XXZ-model in 
Ref. [H]). The later determines the so-called Q-function Q{9, k\ R), whose asymptotic expansions 
at — )• +00 and 9 ^ —00 generate the eigenvalues /2n-i and l2n-i (along with the eigenvalues of 
the nonlocal integrals of motion of Ref. [12] ) , respectively. Remarkable observation of [2l[3] is that 
at the special value /S'^ = | these essentially quantum characteristics can be related to the solution 
of the classical nonlinear differential equation (jl.ip . 

Equation (jl.ip of course is the radial equation for rotationally symmetric solutions of 2D sinh- 
Gordon equation. We will argue that for generic similar relation exists to the classical Modified 
Sinh-Gordon equation (MShG) 

a^a^r/ - e^'' + p(z) p(z) e-2'? = (1.5) 
with the functions p{z) of the form 

p{z) = - . (1.6) 

Here a and s are real, positive parameter^, related to the parameters /3, /u in (jl.2p as follows 

(1-7) 



R 



^7rr(l - /?2 
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but z, z are formal variables, not related to the space-time coordinates (x, t) in (|1.2p . Equation (jl.5p 
with general p{z) is well known in differential geometry (see [13] for review). The same equation 
with polynomial p{z) has emerged lately in different contexts in SUSY gauge theories |14fll8j . and 
these later papers have inspired to large extent the work presented here. 

Obviously, MShG equation in general has no rotational symmetry. Instead, it has discrete 
symmetry 

ITT _ 12L — / \ 

z—f-ei^z, )-e a z. (1.8) 

We will consider solutions of MShG equation ()1.5p which respect this symmetry, continuous at all 
finite nonzero z, z, and grow slower then exponential at z, z — t- 00 (precise conditions are listed in 
Section 2). There is one-parameter set of such solutions, characterized by the behavior r] — )• I log{zz) 
at z, z — )• 0, with real / S [— ^, ^ ] which will turn out to be related to the quasi- momentum as 

/ = 2|fc|-i. (1.9) 

As is well known, the MShG equation is integrable, and the associated linear problem has the 
form |T9] 



D* = 0, D^ = 0, (1.10) 

where D, D are components of sl{2) connectiorH 

D = d, + ^ d,ri - [ ^7+ e" + a" p(z) e''^ ] , (1.11) 
D = 9^7? a3-e-^[cj- e^ + £7+^(2)6-" ] , 

^In fact, for technical reasons in this work we limit our attention to the case a > 1, which corresponds to < i 
in p.2|) . However, our main results remain valid at any positive a. 

2 a , r. ,■ . • ■ 3 A \ + /O 1\ _ /O 

a are the usual Pauli matrices, i.e., ( g I ' ~V0 0/'"" ~\1 
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with the spectral parameter 9. The Q-function of the quantum sine-Gordon model (|1.2|) will 
be related to connection coefficients for certain solution of the linear problem (jl.lOp . Therefore 
our result can be regarded as generalization to ^ 7^ of the relation j20j between the integrable 
structures of CFT [2T] - [23] and spectral characteristics of linear ordinary differential equations [2^ 
125] . Indeed, our derivation follows very closely the analysis in [26]. The novel feature of the massive 
case is that the coefficients in the linear problem are not elementary functions but rather solutions 
of integrable nonlinear partial differential equation (jl.Sp . 

Similar relation exists between certain solution of the MShG equation (jl.Sp . with p{z) of the 
same form (jl.6p but this time with a < — 1, 



a = —b — 1 



^7rr(l + 6^) 

r(-62) 

and the vacuum Q-function of the quantum sinh-Gordon model 



2+26^ 

(1.12) 



{dt^f - {d^y^f - 2n cosh (by^) , (1.13) 



^ 167r 

where again the spatial finite size geometry ()1.3p is assumed. Of course, physical content of the 
quantum sinh-Gordon model is much different from the sine-Gordon model (jl.2p . In particular, 
(jl.lSP has unique vacuum. Correspondingly, the MShG equation (jl.Sp with a < —1 has unique 
solution which is continuous at all finite nonzero z, z. The vacuum Q-function of (|1.13p |27fl29] will 
be related to the linear problem (jl.lOp associated with this unique solution. 

The paper is organized as follows. In Section[2] we discuss the MShG equation with q > 0. 
We define a family of regular solutions, and describe their basic properties. We also discuss the 
associated linear problem (|1.10p . and define the functions Q+{6), Q-{6) as certain connection 
coefficients. In Section[3]we describe how the function Q{9, k) is constructed out of these coefficients, 
and list its basic properties. In particular, we show that it is determined by unique solution 
of complex nonlinear integral equation identical to DDV equation, and thus coincides with the 
vacuum Q-function of the sine-Gordon model. We also establish relation between the classical local 
IM of MShG, and vacuum eigenvalues of the quantum IM of the sine-Gordon model. In SectionH] 
we define the functions Tj{9, k) in terms of the monodromy of the linear problem (|1.10p . and show 
that they coincide with the vacuum T- functions of the model (II. 2p . In particular, at integer values 
of 2a these functions can be determined through the finite system of TBA equations, which were 
previously derived in similar context in |14H18j . The relation is similar to that described in Ref. [30] 
in the massless case. Section[5]is devoted to analysis of the inverse scattering problem in (jl.lOp . In 
particular, we present explicit series-like representation of the solution 77. MShG with a < —1, and 
its relation to the quantum sinh-Gordon model p.l3p . is discussed in Section[6j 



2 MShG equation and linear problem 

Although generally z and z in (II. 5p can be regarded as independent complex variables, in the 
present discussion we usually (but not always) assume them to be complex coordinates on 2D real 
space. Thus, t] in (jl.Sp is assumed to be a function of two real variables, r] = r]{p,(l)), the polar 
coordinates associated with {z,z), 

z = pe"f', z = pe-'^. (2.1) 

We consider special family of solutions of (jl.Sp . parameterized by real / G [ ~ ^' ^] ' defined by the 
following properties: 
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i) Periodicity 

V{P,^+^)=V{P,^)- (2.2) 
In other words, the solutions rj{p, (j)) are single- valued functions on a cone with the apex angle ^, 

C^: </>~0 + f, 0<p<oo. (2.3) 

ii) ?7(p, (p) are real-valued and finite everywhere on the cone C^l, except for the apex /9 = 0. 

ex. 

ill) Large- /3 asymptotic form 

r/(/9, 0) = a log /9 -|- o(l) as p — )• oo . (2.4) 
iv) /) — 7- asymptotic form 

, , f2nogp + 0(l) for \l\<\ , , 
[± logp + C'(log(-logp)) for / = ±i 

Unless specified otherwise, we will describe the cone ()2.3p by the chart, 

Mf^: -q^<4><^. 0<p<oo (2.6) 

with the rays = — ^ and = ^ identified. We assume that solution satisfying conditions i)-iv) 
is unique, and hence respects all symmetries of Eq. (jl.5p . In particular. 



r/(p, (/.) = 7?(p, -0) . (2.7) 
Starting from the asymptotic form (j2.5p one can develop z, z — )• expansion of the form 



ry = / log(zz) + % + X^7fc (^'"' + ^'"') (2-8) 

fc=i 

(1-202 ^ ^ + (l + 2/)2 ^^^^ ' 

where r/o and 7^ are integration constants. It is easy to see that the coefficients in all omitted 
terms in this expansion are uniquely determined once these integration constants are given. On 
the other hand, r/o and 7^ are not new parameters of the solution; they have to be determined from 
consistency of this expansion with the remaining conditions i)-iii)- We will give explicit form of 
the constant 7/0 in Section[3] below (see Eqs. p.l4p . (|3.40p ). 

The expansion (12. 8p remains valid if we regard z and z as independent complex variables. For 
our analysis, the most important message from ()2.8p is that in the "light-cone" limit z — >• (with 
fixed z) 

7? ^ / log(2;z) r/o 7(-2) , (2-9) 
where 7(2;) = X^fc^i 7fc -^^"'^ decays as 

7(z) ~ (2.10) 

at small z. 
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Although the solution rj^p^cj)) is a single- valued functions on the cone (j2.3|) . the connection 
(jl.lip is not. Instead, the linear problem (jl.lOp is invariant with respect to the operation 

0: 4>^^+l, e^d-'^, (2.11) 

involving the shift of the spectral parameter 6. Another easily established symmetry of this linear 
problem involves the operation 

n : e^e-i-K, (2.12) 

which transforms the connection (jl.lip as 

tl[D] = D , ti[D\ = D . (2.13) 

Motivated by these mutually commuting symmetries, we define two solutions ^± = ^±{p, 01^) of 
the linear problem (jl.lOp uniquely specified by their asymptotic behavior 

In writing (|2.14p we have assumed that |/| < ^. In the analysis below we usually adopt this 
limitation, and treat the case |^| = ^ by continuity. Using Eqs. (j2.1ip - (j2.14p . and the fact that at 
real 9 

D{e) = a^ D*{-e)a\ (2.15) 

where the star denotes complex conjugation and = a'^ + a~ , it is straightforward to establish 
the following properties of these solutions: 

• ^±(p, (j) I 6) are entire functions of for arbitrary real (j) and p > 0. 

• 0-invariance: 

*±(p,0 + i^l^-|^)= ^±{p,(t^-^\e + ^) . (2.16) 

• Il-transformation 

*+(p,(/.|6'±i7r) = -e^^^' cj^ *+(p,0|0), (2.17) 
*_(p,(/.|0±i7r) = e^''^' *_(p,0|e) . 

• Normalization condition 

det(*+, *_) = ^— . (2.18) 

cos(7r/) 

Here and below (^+, ^-) stands for the 2x2 matrix with the columns ^+ and 

• For real 9 

^l{p,(t>\9) = a'^^{p,<P\ -6) , (2.19) 

and 

^l{p,0\e) = ^±{p,0\9) . (2.20) 
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Note also that for / — )• zt^: 



lim r cos(7r0 (*+-*-)]= . (2.21) 



- 2 



The above solutions ^± are specified by their p — t- behavior (j2.14p . On the other hand, at 
large p the WKB analysis applies. Assuming that 6 is real, it is straightforward to show that while 
generic solution of (ll.lOp grows exponentially at /) — >• oo, there is a solution which decays in the 
wedge 

- < < . (2.22) 

We denote this decaying solution as H = S{p,(j)\6). It is uniquely specified by the asymptotic 
condition 

S,{p,ct>\e) S-{p,ct>\e) as p-^+oo, (2.23) 
where £- is the shorthand for the decaying exponential 



e 2 \ 

icup exp 

—6 2/ 



^ — cosh.{e + i(a + l)4>) 



(2.24) 



a + 

Since ^± form a basis in the space of solutions of linear problem (jl.lOp . we have linear relation 

H = Q_ {0) *+ + Q+{e) , (2.25) 

where the coefficients Q± {9) (of course independent of the variables p, cp) are functions of the 
spectral parameter 6 as well as the parameter I (the last argument is temporarily omitted in the 
above notations). These coefficients are to be related to the Q-function of the quantum sine-Gordon 
model ([L2|). 

As is well known, the matrix linear problem (ll.lOp can be reduced to second order linear 
differential equations. One can write general solution of (jl.lOp as 

* = f . f \ ] = M i(9,+f^)A , (2.26) 

\ e 2 6 2 [Oz + Ozf]) ip I \ e2 e 2 ip J 
where il) and ^5 solve the equations 

[dl-u{z,z)-e^' p{z)] V = , (2.27) 

[a|-^2(z,z)-e-2%(z)] v; = 0, (2.28) 

with 

u{z, z) = {dzV? - dh , z) = {d-zT^f - dlri . (2.29) 

This form is convenient for making connection to the analysis in Ref. [26] which emerges at 
small s and large 9. Concentrating attention on ()2.27p . consider first the light cone limit z — )• in 
which r] assumes the form (j2.9p . After that, the limit z ~ s — )• , — )• +oo can be taken, with the 
combinations 

X = e^ z , E = s^°' e^ , (2.30) 
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kept finite. Then the term 7(2;) in (|2.9p can be dropped, and (|2.27|) reduces to the Schrodinger 
equation 



+ X 



ij; = E lp 



(2.31) 



In this double limit the solutions ^± assume asymptotic form 



ae(2i+l) 

e 2(a+i) 



/ e^I^ ip+{x) 



{21 + l)Y/cos(7r/) 



\(zz)-2e (9^ + 

t^'+i) / (zz)i e2(^ ^-(2;) 



g 2(ti+l) 

cos(7r/) 



(2.32) 



where ■0± are unique solutions of the above Schrodinger equation (j2.3ip defined for |Z| < | by their 
X — >• behavior 



The equation (j2.25p then reduces to 



ip- — )• X 



-I 



where x is the decaying solution of (I2.3ip . defined by the asymptotic condition 



X — )• 2; 2 exp 



X 



a + l 



+ X 



l-a\ 



as X — )• +00 



(2.33) 



(2.34) 



and 6s = —(a + l) log(s). The coefficients Q^^^^\6) coincide with the Q-functions (the vacuum 
eigenvalues of the operators Q±(^) defined in |22ll23j ) of the left-moving chiral sector of the CFT 
emerging in the massless case /i = in (jl.2p . The proof of this statement is given in Ref. |26] . 
Eq. (|2.25p generalizes (|2.33p to the massive case ^ 7^ 0. 



3 Q-function 

3.1 Basic properties 

The following properties are established by arguments nearly identical to those presented in Ref. |26j 
(For completeness we sketch the derivations in Appendix A) 



• Q±{G) are entire, quasiperiodic functions of 6, 

• At real 6 

Q*±i0) = Q±i0) > Q±{0) = QTi- 



in '\ 
2a J 



(3.1) 



(3.2) 



^Of course the Q-functions of the right-moving sector show up in the opposite double hmit z ^ and z ^ s ^ 

z, -00 in (1^:^ . 
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Q±{9) satisfy the Quantum Wronskian relation 



Q+{e + ^)Q 



2a' 



2a 



2a' 



-2i cos(7rZ) 



(3.3) 



lim Q+{e) = lim Q_(6') . 



(3.4) 



We now define the function Q{9,k) of complex and real k as follows. For k G (— f 2)\^ 



set 



Q{e,k) := { 



Q+iO)\i._ 
Q-{0)\i._ 



:2k- 



'2k- 



for < A; < ^ 
for - i < A; < 



(3.5) 



Due to theproperty (13. 4p this definition extends to /c = 0, /c = it^ by continuity, and then admits 
continuouci periodic extension to all real k by 



Q{9,k) = Q{e,k + l) . 
Thus defined, Q{9,k) is entire function of 9, satisfying 

Q^0 + illMll^k)=e'-''Q{9, k), 



(3.6) 



(3.7) 



Q*i9,k) = Q{9*,k) , Qi-9,k) = Q{9,-k) 



and 



Q{e + ^,k)Q{9 



2a' 



k) - Q{9 - ^, k) Q{9 + ^, -k) = -2i sin(27rA;) . 



(3.8) 



(3.9) 



To fix the function Q{9, k) uniquely, we will need two additional analytic properties. One 
concerns with the asymptotic behavior at |Ke {9)\ — )• oo, another determines the general pattern of 
zeros of this function in the complex 0-plane. Due to the quasi-periodicity (jS.Th one can concentrate 
attention on the strip 



H : \^m{9)\ < . 



Define also 



< Qm{9) < ; 



H. 



■n-jl-i-a) c> 



< 9m6l < 



Then 



For real k and a ^ 7^ 1, 3, 5 

Q - 



r(l + a) 



re^ 2a 



Q e^'^'^Q-^ exp 



4cos(^) 
re 



2a' 



4cos(^; 



£H±, 5Re (9) +00 



eH±, (0) ^ -00 



*Our analysis in subsections l3. 211531 suggests that Q{6,k) is in fact analytic at all real k. 



(3.10) 
(3.11) 



(3.12) 
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Here 



r = B s^^- , B = '^3^^!^^ > (3-13) 



and & is related to the constant tjq in Eq. (j2.8p as follows: 

® = T{l-%) (0 < < ^) • (3.14) 

(3 = ©(/c) is a real function of real variable k, such that 

e{k)e{-k) = i, e{k + i) = e{k). (3.15) 

For any real k, all the zeros of Q{9, k) in the strip H are real, simple, and accumulate towards 
— )■ ±00. Let 



e{e) = i log 



Q[0+'^,k) 



where the branch of the log is fixed by the condition 



(3.16) 



e(^) ^-Wt ^ -2vrA; for Ke(0) ^ +00 and |9m(0)| < ^ . (3.17) 

2cos(2^) 

Then the zeros 9n can be labeled by consecutive integers n = 0, ±1, ±2, . . . , so that On < 
9n+i, and 

e(0„) = TT (2n + 1) . (3.18) 

The asymptotics ()3.12p can be established through straightforward WKB analyses of the linear 
differential equations (jl.lOp . The second property is derived in the next subsection. 



3.2 Zeroes 

Consider first the limit of small s. Due to the relations ()3.6p and (jS.Sp we can assume, without 
loose of generality, that < A; < ^. In this case the problem reduces to Schrodinger equation 
(j2.3ip . and the limiting pattern of zeros of Q{6,k) can be read out from the pattern of zeros of 
Q^^^\9), which is relatively well understood j221l26j. Simple analysis shows that in this case the 
zeros {On} split into two widely separated groups, the "positive" zeros {On}\^=Q, and the "negative" 
zeros {0-n-i}\'^=o- The limiting behavior of the positive zeros follows from Eqs. (l2.30p - (l2.33p . 

^ ^ (3.19) 

where {£'n(0} eigenvalues of Schrodinger operator (j2.3ip associated with the spectral prob- 

lem 

ij; as X ^ (3.20) 

'(/'—)• as X ^ +00 , 
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with I = 2k — ^. If / > — ^ the spectral problem (I2.31|) . ()3.20p corresponds to self-conjugated 
Schrodinger operator, and the associated eigenvalues are non-degenerate and positive. In fact, it is 
possible to show that {£n{l)}^—Q is a family of meromorphic functions of complex I. Furthermore, 
they are analytic in the strip ^e{l) > — | and £n{l) > for real / > — |. At / = — | we have 
£o{-^) = 0, while all higher £n{l) (n = 1, 2...) remain positive. Therefore at s — )• all the 
positive zeros are real, simple and have the limiting behavior (j3.20p for any k G [O, |] , including 
the ends of the segment. 

The limiting behavior of the negative zeros is slightly more delicate, but still can be described 
by equation similar to (|3.19p . 

s2« e-'"'+"r' ^ £:„(-! -/). (3.21) 

Here the argument —1 — / belongs to the segment [ — |, — |] , where 1 — /) are understood as 
the analytic continuations from the domain 3f?e(/) > ^- It is important that all <£'n(— 1 — remain 
strictly positive, except for £q{—1 — I), which turns to zero at / = ^. Thus the limiting behavior of 
the negative zeros is similar to the positive ones, i.e. they are also real and simple. If s is small 
but finite, at A: — )• ^ the single zero 6-i leaves the group of negative zeroes and travels toward 
the positive ones. In the process, 9-i remains real for any < A: < ^, since Q{9, k) is a real 
analytic function. For k = ^ and arbitrary s > we have 9-i ~ 1, and hence it still formally obeys 
Ea. ffsrni) . 

By continuity in s (which we assume) this pattern of well separated groups of positive and 
negative zeros persists for sufficiently small values of this parameter. While exact locations of the 
zeros change with s, no additional zeros can be generated, since this would violate the asymptotics 
form (j3.12p . which is valid for any s. One can define two positive and non-degenerate sets of 
numbers {En{±k)}\'^^Q by 



2ae 



Enik) , n > 

s^" EZl,_i{-k) , n < 



In the domain 0<A;<i, ats— )-0 the sets {£'„(-|-A;)}|^g and {£"„(— fc)} |J!)^o converge to 
\^£n{2k — ^)}|^o ^-^d \^£n{—2k — i)}|^Q, respectively. But at any s we have 

En{±k) ^ {^{2n±2k + l))^ (l-ho(l)) as n ^ +oo , (3.23) 

where the constant B is given by (I3.13p . This follows from the asymptotic (j3.12p . The entire 
function Q(6,k) can be represented as the product over its zeros. Thus, for a > 1 one can simply 
write 



oo 



n=0 



2ka0 -r-r / S e«+l \ / S e a+i , 



where ^{k) is some ^-independent factor. For < a < 1 one has to include Weierstrass prime 
multiplier to make the product convergent. 
Eqs. (l3.6p - (l3.8p imply the relations 

En{k + 1) = En+lik) , Eo{-k - 1) Eo{k) = s^" (3.25) 

and 

(t{k) = -s-2° ^o(^) + 1) , e:(fc) = • (3.26) 
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The overall normalization factor ^ in (|3.'24p and the "reflection S-matrix" (|3.14|) . (|3.15p also can 
be represented as the convergent products 



-1)^2^ lly/Enik)E„i-k) (2n + l))-^ 



(3.27) 



n=0 



and 



re 



,1E 



47r 



4fca 
' Q + 1 



n°° EJ-k) ±kcy 
/ e(Q + l)(n + l) 



n=0 



En{k) 



(3.28) 



where [k] is an integer part of real number k, is the Euler constant, and we use r from (j3.13p . 
The Quantum Wronskian relation ()3.9p shows that at any real zero 9n 



-ie(dn) 



-1 , 



(3.29) 



where e{9) is defined in ()3.16p . Hence, is an odd integer number which, by continuity, cannot 
change when one changes s. Therefore these numbers can be extracted from the limit s — )■ (see 
e.g. Appendix A in Ref. [31j). By this argument, Eq. (j3.18p is valid for any finite s. 

Since t{9) = + fc is real valued, monotonically increasing function of the real variable 0, one 
can introduce the inverse function 6{t), which has the properties 

e{T) = -9{-t) (9m(r)=0) 

(27rT \ 
j as r — +oo . 

Then, 



(3.30) 



Eoik) 



„2a 



exp 



^,9{2k + l) 



(3.31) 



and 



En{k) = Eo{k + n) , En{-k) 



Eo{k-n-l) 



(n = 0, 1...) 



3.3 Destri-De Vega equation 



The quasiperiodic entire function Q{9, k) is completely determined by its zeros 9n in the strip (jS.lOp . 
and the asymptotic condition (|3.12p . On the other hand, the positions of the zeros 9n are restricted 
by the equation (|3.18p . Mathematically, the problem of reconstructing the function Q{9,k) from 
this data has emerged long ago in the context of the analytic Bethe ansatz |32H35j . For the sine- 
Gordon model (jl.2p . the problem was solved by Destri and De Vega [9], who have reduced it to a 
single complex integral equation, the DDV equation. 

Starting with the equations ()3.24p , ()3.18p , (|3.12p , and following the steps outlined in [9HTT] , one 
can derive the equatioio 

2 / d^' 9m [ log ( 1 + e-''(^'-'°) ) 1 ^(^ - 9') = r sinh(0) - 2TTk - e{9) , (3.32) 

J — oo 



^To be precise, the derivation in [9lll0j goes through the analysis of discretized version of the sine-Gordon system, 
where the analytic properties of the counting function are somewhat different, with the continuous limit taken at the 
end. Only minor adjustment of the DDV arguments are needed in the present context, see [22) . 
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where 



G{e) 



du e'^^ sinh 



' 7n^(l— a) ^ 
. 2a y 



27T 2 cosh (^) sinh(f^) 



(3.33) 



and r is given by (|3.13p . 

Equation (|3.32p coincides exactly with the DDV equation for the sine-Gordon model (jl.2p if 
one sets 



a = /3 



1, 



r = MR, 



where |36 



M 



2jP 



vr/ir(l- /3^) 

r(/32 



-, 



2-20^ 



(3.34) 



(3.35) 



is the sine-Gordon soliton mass, and identifies k with the sine-Gordon quasi-momentum. We 
conclude that the function Q(6,k), defined in terms of the coefficients in (j2.25p as in Eq. (j3.5p . 
coincides with the Q-function |32ll33l[37l - l39j of the quantum sine-Gordon model. 

Few remarks about the equation (j3.32p are worth making. Note that at a = 1 the kernel G{9) 
vanishes at all 6. In this case the solution of the DDV equation takes simple form 



e(6')|a=i = vrs^ sinh((9) - 27rfc . 
One can use p.36p to obtain explicit form for EQ^k) (I3.3ip in this case, 

Eo{k)\a=i = 2k + l + y^{2k + lf + s^ . 



(3.36) 



(3.37) 



For general a > 1, Eqs. (|3.16p and (|3.12p allow one to reconstruct the Q-function. In particular for 
^m{9) = ttis middle of the strip if+ (13. lip ) Q{6,k) is expressed in terms of the 

solution of DDV equation as follows: 



log + A.) = ^^ifl + ivr^ + 1 log(6) 



-2i 



d^' 9m 



(3.38) 



2a ^ 



log ( 1 e' 



-ie(e'-iO) 



F{e' -e-io) 



Here, 



m 



-oo 2vr 4cosh(f ) sinh(^i^) 



and 



r-oo 

log((3) = a I — 9m 

1 TT 



log f 1 + e--(^-iO) 



(3.39) 



(3.40) 



Once €{6) is determined at real values of 9, Eq. ()3.38p can be modified to yield the Q-function in the 
whole strip H^, and then by (j3.ip in the whole complex plane. For < a < 1 Eq. (j3.38p requires 
minor modification. 
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3.4 Large-6' asymptotic expansion and local IM 

Using (j3.38p it is straightforward to obtain the large-^ asymptotic expansion of the Q-function. For 
I $5771 (0)| < '^^'^^^ we have 



logQ(0 + M|±i)|fc) 



r e 



-1 e 



-(2n-l)e 



4cos(^) 

oo 

+ E ®" 



— + mk + \ log(6) 



(3.41) 



n=l 



n=l 



as SRe (^) — )• +00, and 



T(a+1) 



2a 



r e 



4cos(^) 



+ i^A; - I log(©) 



(3.42) 



n=l 



,(2n-l)e + Y,^n 
n=l 



2ane 



as (0) — — oo. Here we use the notations, 

'^2n-l = 



rSn,i ^ (-1)"+! 



4cos(^) ' sin(^(^: 



d9 
vr 



3(2n-i)(e-io) (1 ^ e-i.(e-io)) 



and 



'^2n-l 



cos(7ran) 



-oo 



(3.43) 



4cos(^) sin( 



n+l /.oo 



(-1) 



« (-1) 

cosfvran 



n roo 

oo 



2a„(e-i0) i^g (1 ^ g-.«?-iO)) 



(3.44) 



On the other hand, the asymptotic expansions ()3.4ip and ()3.42p can be obtained directly from 
the WKB expansion for the hnear differential equations (ll.lOp . To simplify calculations, it 
convenient to trade the variable z to 



IS 



w = dz \fp{z) , 



(3.45) 



and similarly for w. As is well known, this transformation brings MShG equation (jl.Sp to the 
conventional Sinh-Gordon (ShG) equation 



for 



with 



dn,d^n - e^^ + e'^'' = 



'n = v-j log (pp) , 



(3.46) 



(3.47) 



P = P{z) 



P = P{z) . 



(3.48) 
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In what follows, we will assume that w = w* , and choose the branch of \/p{z) in (|3.45|) which is 
positive at the upper edge of the branch cut in Fig(TJ 

Using the standard technique |19] for the large-^ expansion, one can find explicit expressions 
for the coefficients CJ2n+ii and 312^+11 as the functionals of 7). The coefficients 32n+i and 
^2n+i appear as local functionals, i.e. the integrals of local densities, 



3'2n-l 



2(2n-l)sin(^i^) Jc^ L 



dw P2n + <^W i?2n-2 



and 



J2n-1 



2(2n - l)sin( 



7r(2n-l)^ 



dw P2n + dw i?2n-2 



(3.49) 



(3.50) 



The integration contour here is a w-image of the contour Cz in the z-plane shown on Fig{Tl 
while Cu) = C^. The functions {P2n, R2n-2) and {P2n, R2n-2) = (-Pirn -^2n-2) conventional 



© 




Figure 1: The chart A4^z'\ and the integration contour Cz in Eq. (|3.58p . 
tensor densities of the local IM for the ShG equation (j3.46p , satisfying the continuity equations 

dwP2n = dwR2n-2 ■ (3.51) 

They can be obtained in explicit form as follows. Let 

u = {d.ujfif - dlf] , V = [du,T)f + dl^fj . (3.52) 

Then 

An = Un[u] , ^2n-2 = e~^^ C/„-l[i)] - 5„,l , 

where C/n[u] are homogeneous (grade(ii) = 2, grade(9) = 1, grade(?7n) = 2n) differential polyno- 
mials in u of the degree n (known as the Gel'fand-Dikii polynomials |40j). 

[/„[«]= A"- 1. (3.53) 

Here 

A = -i 52 + u- I n', (3.54) 
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and prime stands for the derivative. Thus, 



Uo[u] 
Ui[u] 
U2[u] 

Un[u] 



1 , 

in , 

3-2 1 '■// 



(3.55) 



16 " 32 / 16 ^ 32 " 



r(n+|) 



J2 

+ . . . 



where the last hne shows overah normahzation of the polynomials. 

Of course, it is straightforward to rewrite the integrals (j3.49p . (|3.50p in terms of the original 
variables z, z, in which the integration contour is just Cz in Fig{TJ We have 



dw = -s/ p{z) dz 



u = p 

where p = p{z), u = {8^1])'^ — d'^r], so that 



u + 



App" — bp' 
IQp^ 



/2 



P2n[u] =P"" P2n[u] + ... 



For example 

3i sin(^) 



7r(2a - 1) 1 
12r(a + l) ^2 Jc^ 



dz — ^ + dz y/p (\/pp e - l) 



(3.56) 



(3.57) 



(3.58) 



Here r is given by (j3.13p . and first term in the r.h.s. is the evaluated form of the integral ^ dz 
p-| [App" - 5p'2). 



3.5 Local IM in quantum sine-Gordon model 

As was mentioned in Introduction, the quantum sine-Gordon model has infinitely many local inte- 
grals of motion, 

l2n.-l 
l2n-l 

where . . . stand for the terms involving higher derivatives of if, as well as the terms proportional 
to powers of /i. The displayed terms fix the normalization of these operators. We will denote 
-^2n-i = hn-iik) and /2n-i = hn-iik) the A:-vacuum eigenvalues of the operators ()3.59p and 
p.60p . respectively. In the CFT limit M = (i.e. at = in ()1.2p ) these functions become 
polynomials in k of the degree nlfl, and the normalization in (I3.59P is such that at M = we have 

/2n-l(fc)L,=o = ^2n-l(fc)L,=0= (f)'""' [ (^Z^)'" + • • • ] ■ (3-61) 

The expansions (|3.4ip . (|3.42p are in agreement with the expected asymptotic behavior of the 
Q-function of the quantum sine-Gordon model |22y28j. with the coefficients 32n-i, ^2n-i and ©„, 

® In this limit (|1.2|) acquires continuous symmetry with respect to any shifts of the field (p, and the analytic 
continuation in k is no longer periodic. 



-4n 



-4n 



R 



dx 



2^ 
^ dx 

2^ 



(3.59) 
(3.60) 



16 



&n being (up to normalization) the fe-vacuum eigenvalues of the local and non-local [12] integrals 
of motion. In particular, we have 



'^2n-l — Cn k 



2n-l ; 



^2n-l — Cn k 



2n-l 1 



(3.62) 



where Cn are constants, independent of k and R. Their exact values can be found by comparing 
the s — )• limit of 32n-i to (I3.61[) . 



Cn 



_ 2n-l ^ p^ (2n-l)(a+l) - 



Here 



m 



a+lJ 



2M sin (^) 



m 
87^ 



r(^)r(-^) 



l-2n 



(3.63) 



(3.64) 



At a > 1, this quantity coincides with the mass of the lowest soliton-antisoliton bound state of the 
sine-Gordon model. 



4 T-functions 

4.1 Definition and relation with Q± 

Let us consider the action of the symmetry transformation Cl, Eq. (j2.1ip . on the solution (j2.23p . As 
follows from (j2.25p and analyticity of ^± and Q±, the solution H is entire function of for any 
real p and (f). Therefore analytic continuation of H in can be used to specify another solution of 
the linear problem (|l.lUp . 

Hi(p, I 0) = 0[H](p, I ^) = H(p, <A + ^ I - f ) . (4.1) 

It is easy to see that for real 9 and < 2(a+i) '^^ grows at large p as 



Siip,^\e)^£+{p,^\0) , (4.2) 



where 




exp 



+ — cosh(6l + i(a + 1)0) 

a + 1 ' 



(4.3) 



Since 



det (H, Hi) = det (5_, 5+) = -2i , (4.4) 



the pair H and Hi forms a basis in the space of solution of the linear problem (ll.lOp . 

Furthermore, applying the symmetry transformation fi" with n E N, one can generate an 
infinite series of solutions 

j]"[H] ^H(p,0+^|0-if^) . (4.5) 

Of course, each of these solutions is a linear combination of the basic solutions H and Hi = 17[H]. 
Using Eq. (l2.25]) . it is straightforward to show that 

0"[H] =-r^(e-i^) H(p,</.|e)+T^(0-^) Hi(p,<^|0) , (4.6) 
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where 



Tm - 



2 cos(7r/) 



+ 



i7r(2i + l)^ 
2a 



Q-{0 



i7r{2j + l)^ 
2a I 



i^(2j+l) ^ 
2a y 



(4.7) 



i7r(2i + l)\ 
2a I 



Note that 



r_i(0) = O, ro(e) = l, (4.8) 



where the second identity follows from the Quantum Wronskian relation (j3.3p . 

Tj(0) can be interpreted as Stokes coefficients defining the large p-behavior of the basic solution 
H. Indeed, Eq. (|4.6p shows that for real 0, and (j) in the domain 

2{a+l) <- 9 2(a+l) ' n - i, Z, . . . [^.y) 

the asymptotic of H can be described as follows, 

s{p, <p\e)^ (-i)[^] {e + ) 5-.(/9, ^\e)+ (4.io) 



rn-li 

2 

(-l)[t] Tn^(6/ + i^) £^(/5,0|0) with a = (-l)" as p ^ oo 



Here are given by Eqs. (j2.24p and (j4.3p . Of course, for given integer n only one term in ()4.10p 
is relevant whereas another term should be ignored. 



4.2 Fusion relations and F-system 

The coefficients Tj{6) in ()4.5p have useful interpretation in terms of the quantum sine-Gordon 
model. Tj{6) are the /c-vacuum eigenvalues of the commuting "transfer-matrices" Tj{9), the traces 
(over the 2j -|- 1 dimensional auxiliary spaces) of the quantum monodromy matrices |33p34| (see 
also Ref. [41j, and Appendix B). The analytic properties of Tj{6), which follow from (j4.7p . are 
in agreement with expected analyticity of the Tj(0)-operators in sine-Gordon model [H]. The 
identities 

T, {6) T, {6 + M|±i) ) = r,._ . {e + M|±H) ) + T,.^ , (0 + (4.11) 

(another simple consequence of (j4.7p ) coincide with well known fusion relations for the sine-Gordon 
transfer- matrices . 

The fusion relations can be taken as the starting point in purely algebraic derivation of the 
TBA equations. In the sine-Gordon model with generic the TEA leads to an infinite system of 
coupled integral equations, the Takahashi-Suzuki system [l2], which is somewhat difficult to deal 
with. However, at rational values of it truncates to a finite system. Thus, at 

a = iV- 1 = 2, 3, ... or a = f - 1 = i, |, . . . (4.12) 

Eq. (|4.7p dictates additional relation 

TN{e) = 2 cos{2TTk) +TN_^{e) , (4.13) 
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which closes the fusion relations (|4.11|) within a finite number of functions Tj{0), j = ^, 1, |, . . . y — 
1 (This truncation is discussed in Ref. ^41j ) . The standard way to derive the associated finite system 
of TBA equations is to introduce the functions 

Y,i9) = r,._iWr,.+|W (j = i,...,f-i) , 

Yo{0) = 0, (4.14) 

Y{e) = TM^,{d). 

As follows from the fusion relations and Eq. (l4.13p . Yj{9) satisfy closed system of functional equations 
(the so-called y-system) 0: 



+ f ) - 1 ) = ( 1 + (^) ) ( 1 + (^) ) ' J 



1 Af-3 
2' • • • ' 2 ' 



-1 + i) -1 - i) = ( 1 + >^iv_3 (^) ) ( 1 + e^-'^ Y{e) ) ( 1 + e-^-'^ Yi0) ) , 

^(^ + 5l)^(^-il)=l + >^f-iW- (4-15) 

This truncated Y-system coincides with the functional form of the TBA equation of Dat type [8] and 
can be transformed to a set of the integral equations. This gives an alternative way to reconstruct 
all the functions Tj{0) and Q±{9) in the case integer 2a ()4.12p . 

4.3 Basic properties of Ti 

2 

In what follows, the function T{9, k) = Ti{9) plays special role. For future references, let us sum- 
marize here its general properties. All the statements listed below are straightforward consequences 
of the definition the definition (14 .TP and the properties of Q-function. 

• For real k, T(9, k) is an entire function of 9, even and periodic in this variable, 

Ti-9,k)=Ti9, k), 

T(0 + ilMl), A;) =r(e,A:) . (4.16) 

• For real k, T(9, k) is a real analytic function of 9, 

T*{9, k) = T{9*, k) . (4.17) 

• It is even periodic function of k, 

T{9, k) =T{9,-k) , T{9, k + 1) = T{9, k) . (4.18) 

• It satisfies the Baxter's T — Q equation 

Ti9, k) Q{9, k) = Q{9+'^,k)+Q{9-'^,k) . (4.19) 

• For a > 1 and Ke(6') — )• ±00 



T{9, k) exp 



4 f 

— cosh(0) 



cos V 2a 



for 1 9m(e) I < . (4.20) 



Note that the TBA equations (|4.15|l for I = 0, i.e. e^"'* = i in (|4.15|l . correspond to the solution of MShG which 
remains finite at the apex of the cone Cil (|2.3[) . This case is of special interest for the problem considered in [161117) . 



19 



Here and bellow we use the notation, 



" s'^" = ^ , (4.21) 



(a + i)r(-^)r(^) ^''^^ 

where m is the same as in p.64p . 



5 Inverse scattering problem for MShG 

5.1 Gel'fand-Levitan-Marchenko equation 

As follows from the asymptotic formula (j2.23p . the composition of symmetry transformations (|2.11|) 
and (|2.12p 11 o J^, acts irreducibly on the solution H, 

Uon[s] :=s{p,<p+'^\9-'-^^) = -ia'S{p,c^\e) . (5.1) 
Combining this equation with (14. 6p one obtains 

( -ia' r 3(P, I ^ + ivrn ) = -T^ {6 - '-^^ ) S{p, <t>\9)- (5.2) 
ir^(e_i™)^3 H(p,0|^ + i^) (n=l, 2...). 



For n = 2, with Tq = 1 and the periodicity condition for Ti, Eq. (|4.16p . taken into account, Eq. ([^ 

2 

becomes a simple difference equation, 

S{p,cl)\e + m)-S{p,(l)\e-i7r) = [Ti{0) S{p,(j)\e) . (5.3) 



Our next goal is to transform (15. 3p into an integral equation defining H. Let us note that the 
coordinates (p,(j)) on Cn, (12. Sp . appear in (15. 3p as parameters. For our analysis here, it will be 
convenient to use slightly different coordinates on Cil. We redefine the angular coordinate (j) by 

a 

shifting it by half the period, 

+ (5.4) 
and use, instead of A4^^\ the chart A4^z^^ defined by the same conditions as in (|2.6p but in terms of 



the shifted angle (j). Correspondingly, z = pe^''' will now stand for the rotated complex coordinate, 

z — > e2^ z . (5.5) 

After this rotation we have p{z) = — (z-^" + s^"). The minus sign here does not affect the form of 
p.Sp . and in this and subsequent sections we set 

p{z) = + . (5.6) 

In the rotated coordinate the zero of p{z) appears at the point (/j, (/>) = (s, ^) ~ (•^'"iz^)' 
shown in Fig. 2a. Actually, we will need the coordinate w on related to the new z as in (j3.42p 
(it differs from w in ()3.45p essentially by a phase factor). To be precise, we set 



yj = w{z) = ^^-^+ I dC VC^+^-C ■ (5.7) 
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2f 



The image of the chart A^i"*"^ in the tw-plane is shown in Fig. 2b. In M^w'^ the apex of the 

cone has the coordinates (^0,1^0) with 

Wo = w{0) 

while 

;(se'^^) = -f cot (^) ± i 



'sin(S) ' 



(5.8) 



(5.9) 



where notation ()4.2ip is used. In the remaining part of this paper we discuss in terms of these 
redefined coordinates z (or w), unless stated otherwise. 
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(a) 




Figure 2: (a) The chart Ai^2^\ Two boundary rays (solid lines) are identified. The dots indicate 
position of the zero of ^(2:). (b) The image A^L*"^ of Ai^^^ under the map ()5.7p . The segment A'B' 
is identified with A'B', and the boundary line from B' to infinity is identified with the line from 
B' to infinity. Specific regions in this chart are discussed in the text. 



In the new variables {w, w) the solution H(p, (p\0) of the linear problem ()1.10p can be written 



as 



[w, w 



H(p,0|e)=exp[^ (i^ + log(|)) 
Here H('u;, w\6) \s the solution of the linear problem 

dy, + \ dyji] a'^ - [o-+ e'^ + cr"e" 



i7r(a+l) 
2a 



9,7 



(5.10) 



(5.11) 



associated with ShG equation p.46p . satisfying (for real 9) the asymptotic condition, 



-1 



exp ( 



w e 



w e 



') 



as w 



(5.12) 



Note that S{w,iv \ 9) is nothing but conventional Jost solution for (|5.1ip (see, e.g., [TO] for details). 
The main advantage of the coordinates [w, w) is relatively simple form of the large-0 asymptotic 
behavior of H. Simple analysis of the linear problem ()5.1ip shows that 



S{w,w\ 9) exp (we^ + tfe ^)— ^ e'^^ e^ ''2 



1 



as 



9 ^±00. 



(5.13) 
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where lo± = u}±{w, w) are local solutions of the Laplace equation d^dyjUj^ = (in fact, u:± 
are piecewise constants, see Eq. (j5.22p below). Combining the difference equation (j5.3p with the 
asymptotic behavior ()5.13p and the analyticity condition for H, it is straightforward to show that 



l{w^ w I 6) 



Kj^Xj^{w, w I 9) 



exp (we^ + we ^) , 



-K^X-{w, w I 9) 

where k± = K±{w,iju) do not depend on 9, while X± are solutions of linear integral equations 

hoc 



(5.14) 



X±{w,w\9) = 1± / 



with the kernel 



tsinh{^) D{w,w\ 9') X±{w,w\ 9') 



D{w,w\9) = Ti {9 + ^^^) exp {-2we^ -2w e"^ ) 



With Eos. dSl^ . (imil and 

det(s{w,w\9 + f), S{w,w\9-f)) = 2 



one obtains 



= d 2 



1-dl 



where 



d± 



+ 00 



d9 



D{w,w\9) X±{w,w\9) 



and 



X+ [9 + f )X_ - f ) + X+ [9 - f )X_ (0 + f ) = 1 - d+d 



(5.15) 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



(5.20) 



In these equations and below we omit the argument {w,w) in X± and d±^ n±. Comparing the 
9 — >• iboo limits of ()5.15p with ()5.13p , one can express the solution 17 of ShG equation ()3.46p itself, 
as well as uj± , in terms of d± , 



(l + d+)(l + d_) 

(l-d+)(l-d_) 



and 



1 ± 



d^ — d_ 



d 



Note that the conjugation condition 

a*{w, w\9) = —a^ 'S{w, w \ — 9*) 

implies the relations 

X_{w,w\9) =Xl{w,w \ -9*) , 



K- = K, 



+ ' 



(5.21) 



(5.22) 



(5.23) 



(5.24) 
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and hence 



Eq. (|5.16p can be rewritten in the form 



D{w, w\9)= rf (0 + ilMll) exp ( - 2ws e' - 2ws e"^ ) 



(5.25) 



(5.26) 



where (see Eq. (j5.9 



Ws{z) 



2a 



+ S 



2a 



W[Z 



se ^ 



+ f tan(^)-i 



and 



T["'^\e) = Ti{e) exp 

5 2 



4f 



cos 



cosh(0) 



2a' 



(5.27) 



(5.28) 



Advantage of using the function T^^°^^ (9) is that it is weh defined at any a > 0, including the points 
= 1, 3, . . . , as opposed to both Q±{9) and T{9), which at these points are defined only modulo 
overall factor exp( const 0). Eq. (l5.26p makes it explicit that the kernel of integral equations (15.15P 
is well defined at any a > 0, including the integer points. For example, in the case a = 1, it follows 
from Eqs.dSlS]), (I4l9]) and (lOHD that 



exp 
log 



_^ 27rcosh(6l - t) 

-j^ _j_ g— TTS'^ cosh(t)+27riA; ^ -j- g"'"'*^ cosh(t) — 27rifc ' 



(5.29) 



For k = ^ this function, interpreted as a Stokes coefficient, was found in Ref. [n]. This result was 
also used in Ref. [lejl. 

Eqs.(l526D-(l52H]) and ([TO]) imply that for 



5Re {w) > - f tan 



(5.30) 



and real 9 the kernel \D{w,w\9)\ is bounded by exp ( — C cosh(^)) with positive constant C. 
In this domain of w the integral equations (|5.15p have unique solutions obtainable by iterations. 
With appropriate deformation of the integration contours in (j5.15p one can extend the domain of 
applicability of the iterative solution to the region /' C M.^^"^ shown in Fig. 2b. The most efficient 
way to solve ()5.15p at the boundary of /', i.e. on the segment 5Re {w) = —f tan (^), {w)\ < f, 
is based on the integral transformation in A = e^ generated by the kernel exp ( i,^ (A — A~^) ) . The 
transformation brings (j5.15p into the form of Gel'fand-Levitan-Marchenko equation (see e.g. [19] for 
details). Important problem of reconstruction of the Jost solution S{w, w\9) in the whole domain 
Aii^^ is beyond the scope of this paper. 



^In the notations of Ref. [T4l[T6]: 71(C) = T("'°s)(6I) with Ce" 
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5.2 Solution of MShG equation 



Since fj decays at — )• oo, at large \w\ it approaches certain solution fji of the linearized equation 
dwdw fji — Afji = 0. The form of this solution can be read out from Eq. ()5.2ip . 



i7r(a+l) 
2a 



) exp(-2'u;e^ -2u;e"^) 



(5.31) 



Furthermore, iterations of (|5.15|) produce (through (|5.19|) . (|5.2U|) . and (|5.21|) ) systematic large-|t(;| 
expansion of r). In fact, one can guess the form of this expansion without explicit calculations. As 
is known [33], [3], the series 



2 



n=l 



mn-1 



2n - 1 



2n-l 

n 



d9j r{dj) exp {-2we'^ -2w e"^ ) 
4^ cosh C-^^^) 



(5.32) 



with = 01, and arbitrary function T{6), provides formal solution of the ShG equation (j3.46p . 
Since the asymptotic form ()5.2ip fixes the solution uniquely, we conclude that the solution fj we 
are interested in is given (in certain domain of w specified below) by the series (I5.32p with 



(5.33) 



This representation is very useful since at sufficiently large w the series ()5.32p converges fast. In 
view of Eq. (j4.20p . the multiple integrals in (15.32P converge only when w belongs to the domain 
(j5.30p . With deformation of the integration contours, the convergence domain can be extended to 



/' C M^j~^ shown in FigiSb. Thus, the solution of MShG ([LSI) in the domain / C Ml'^' (FigiSJi) can 
be written in the form rj = fj + j log{pp), where fj and p are given by (I5.32p and (15. 6p . respectively. 

It would be interesting to find similarly explicit expression for r/ in the remaining part of 

As {wl — )• oo, the integral (j5.3ip can be evaluated by the saddle-point method. As the result 
one derives the large-/? asymptotic of the solution of MShG (|1.5p . (11. 6p . 



(+) 



\ log ( - 2 (sp)2- cos(2a0) + ) 



(5.34) 



+ Ti(i(a + 1)</.) 



V2 



ITT 



1 + 0{t~^)) as T 



a+1 



+00 . 



In writing this equation we use the original polar coordinates (/?, (p) on the chart Ai^^^ (j2.6p . and 
assume for simplicity that a > 1. Eq. (|5.34p shows that the T- function Ti (9) determines the angular 

2 

dependence of the sub-leading large-p asymptotic of the MShG solution. At s = and finite 9 the 
function Ti{9) becomes a constant. 



nWL=o = 2cos(^), 



(5.35) 



and hence (j5.34p generalizes the well known asymptotic formula for the Painleve III transcendent 

mi. 
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6 MShG with a < —1 and quantum sinh-Gordon model 



So far we were concentrating attention on MShG with a > 0. In that case we had a freedom to 
adjust the asymptotic behavior of r/ at z — )• as in (j2.5p . with the free parameter I. If a < — 1, 
the situation is different: asymptotic form of a regular solution is fixed at both z — )• oo and z — t- 0. 
In fact, the roles of z = and z = oo can be interchanged by certain "duality" transformation. 
Conformal transformation 

a+l 

(6.1) 



a + l 



(with suitable shift of the field rj) brings the MShG equation to the original form, but with the 
parameters a, s replaced by the "dual" values a, s related to the original ones as 

(a + l)(a + l) = 1, s-+i = _i!li. (6.2) 

a + l 

Note that in terms of the parameters b and fi, Eq. ()1.12p . these relations faithfully reproduce the 
6 — )• duality symmetry of the quantum sinh-Gordon model |44) . Below we argue that such 
solution (more precisely, certain connection coefficient for solutions of the linear problem (ll.lOp ) is 
related to the Q-function of the sinh-Gordon model ()1.13p . 

In this discussion we will use the chart A^i"*"^ on Ci, in which p{z) has the form (15. 6p . and 
associated polar coordinates z = pexp{i(j)). We also use the notation 

u = —a . (6.3) 

We assume that MShG equation has solution ri(p, (p) which is real and continuous on Ce. except for 
the apex p = 0, and has the following properties (compare to the properties i) — iv) in Section[2]) 
i') 

7?(p,(^ + f) =r?(p,</.). (6.4) 
ii') r](p, 0) are real- valued and finite everywhere on the cone Ci, except for the apex p = 0. 

a 

iii') 

r/(/9, (/))—)• as /9 — )• oo . (6.5) 



IV 



r]{p,(j)) -V \og{p) + 0(1) as p ^ . (6.6) 



The relevant linear problem (jl.lOp now is somewhat simpler then that previously discussed. At 
> 1 we introduce two solutions of the linear problem, defined by the asymptotic conditions at 
p — )• oo and p — )• 0, 



(P, </'l^) 



1 



exp 



2s p cosh(0 + i(i 



as /? — )• +00 



(6.7) 
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and 

3_(p, <p\e) 

We then define 



e 2 

ivcf} 



exp 



cosh(0 - - 



as /9 — )• . 



Q(0) = i det (H+, H_) . 



(6.8) 



(6.9) 



By arguments parallel to the analysis in |44p45j . it is possible to establish the following properties 
of this connection coefficient: 



Q{9) is an entire function of 9 with the symmetries 

Q{o) = Q{-e), Q*{e) = Q{e*). 

Q{6) satisfies the Quantum Wronskian relation: 



2i/ . 



1 . 



(6.10) 



(6.11) 



As the function of the complex 9, Q{6) is free of zeroes in the strip |$5m0| < ^ + e for some 
finite e > 0. 



For |9m0| < f + e and Ke6' +oo 



Q{e) = exp 



2 f 



wher^ 



sin(^) 



3 

7r2 1/ 



^^" + 0{e- 



Let us introduce the function e{9) through the relation 

Q(^ + f)Q(0-if) = l + e-=(^) 



(6.12) 



(6.13) 



(6.14) 



With the analytic properties listed above, it is straightforward to transform the difference equation 
(|6.1ip into integral equation for e{6) (see Ref. [27]). 



with the kernel 



Then 



/oo jnl 
— ^{9 - 9') log (1 + e-^(^') ) 

4sin(^) cosh(0) 



m 



log Q{9) 



4f 



cosh(26') -cos(^) ■ 

d^' log (1 + e-^(''')) 



sm - 



— cosh(6') + 



27r cosh(6' - 6^') 



(6.15) 



(6.16) 



(6.17) 



®Note the similarity in the definitions of f for v = —a > 1 and for v — — q < (|4.21l) . It is defined to be positive 
for any \a\ > 1. 
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As was argued in Ref. [25], Eg . (16.171) gives the Q-function for the quantum sinh-Gordon model in 
the finite-size geometry (jl.Sp . provided f = with m interpreted as the mass of the sinh-Gordon 
particle, and u related to the sinh-Gordon coupling constant b as in ()1.12p . i.e. 



i^ = b-^ + 1 . 

It is instructive to review the above statements in terms of the coordinate 



w = w{z) = dz \l + 5-21^ 



(6.18) 



(6.19) 



which brings the MShG equation to the form of the conventional ShG equation ()3.46p for f] 
V ~ \ log(pp). We fix the integration constant in ()6.19p in such a way that 



( se=''^2^ ) = ±ii 



Explicitly, 



w{z) 



cot (^) + z ,F^{- 1, -± 1 - ^, -(f )-'^ ) . 



(6.20) 



(6.21) 



The chart A^i^^ and its u;-image are shown in FiglS^ and FiglSb. The asymptotic conditions 



<(+) 



© 




(a) 




Figure 3: The chart TWi^^ (a), and its image ^A^' under the map ()6.21[) (b) in the case v = —a > 0. 

iii') and iv') above simply mean that fj decays at \w\ — )■ oo in both regions /' and //' in Fig. 3b. 
Eq. ()6.9p can be equivalently written as 



<(+) 



Q{9) = \ exp 



4r 



sm(^) 



^ cosh(0) 



det ( H+, H_ 



(6.22) 



where H-t are conventional Jost solutions |19] for the linear problem (15. lip satisfying the asymptotic 
conditions 



S+(w, t(; I 0) — )• exp ( — w e^ — w e ^) 



1 



as 'Re{w) — )• -|-oo , 



(6.23) 
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and 



H_(tt;,t(; I 0) — )• exp (we +we ) i-^ 



as Ke(it;) 



-oo . 



(6.24) 



In this picture the duahty ()6.ip is quite evident. The rotation w — t- —w interchanges domains 
/' and //' in FigjSh , which is equivalent to the change of the parameters 



r ^ r 



(6.25) 



identical to ()6.2p (Note that (I6.13P is invariant under the transformation (j6.2p ). It is easy to check 
that under this duality 



{w,w \ 9) ^ H=p(-w, -w I 9) 



(6.26) 



and hence the function Q{9) does not change when its parameters are transformed as in (j6.25p . 
Thus, it is invariant with respect to the 6 — )• duality, as the sinh-Gordon Q-function should be. 

The asymptotic expansions at — )• +oo and 9 — )• — oo of Q{9) generate the vacuum eigenvalues 
of the local integral of motions in the quantum sinh-Gordon model. Namely [28J, 



logQ~-Coe^-^C„/2„-ie 



-{2n-l)9 



as 



+00 , 



(6.27) 



n=l 

oo 



log Q ~ —Co e 



J2Cnhn-ie^'-'^' as 



-oo 



n=l 



where /2n-i and l2n-i are the vacuum eigenvalues of the local integral of motions l2n-i and l2n-i 
(n = 1, 2, . . .) normalized as in (I3.59p . The constant C„ are given by (I3.63P with a = —i', where 
m is interpreted as the mass of the sinh-Gordon particle. From Eq. (j6.17p we have 



Cn l2n—l 
Cn l2n-l 



Co 5n.i + (-1)" ^ e(2-i)^ log (l + e-(^) ) , 
Co 6n,i + i-ir f ^ e-(2"-i)^ log (l + e-W ) . 



(6.28) 



On the other hand, straightforward WKB analysis allows one to express (I6.28P in terms of the 
classical conserved charges for MShG equation. 



Cn hn-i 

Cn l2n-l 



i-iy 



n foo _ 



2n- 1 

(-1, 
2n- 1 



oo 

n /"OO r 



dw P2n + dtf -R2n-2 
dw P2n + dw R2n-2 



(6.29) 



Further analysis is similar to the one presented in Section[5] for a > 0. The novel property is 
existence of two Jost solutions. For this reason one has to introduce two T- functions 7±, 



T+{9) Q{9) = Q{9+'-^) + Q{9-f) , 
T^{9)Q{9) = Q(0 + ii^) + Q(0-Mf^pi)) . 

Note that the duality transformation ()6.25p interchanges the T-functions, 

T±{9)^T^{9) . 



(6.30) 



(6.31) 



28 



The ShG solution i) still can be written as the series (j5.32p . but the choice of the function T^O) is 
different for the two parts I' and II' of the chart A4[^^ in Fig. 3b. Thus, for w ^ I' we have 



V 



E 

n=l 



2n - 1 



2n-l 

n 

i=i 



dOj T+{9j) exp ( - 2w+ e^ - 2w+ e"^ ) 
4^ cosh C-^^i^) 



where 

w+ = f cot (^) + w . 
For w £ II' one has to replace T+{9) by T-{0), and by 

— f tan 



W- 



w . 



(6.32) 



(6.33) 



(6.34) 



Since the union /' U //' covers the whole chart combination of these two representations 

provide the solution i] on the whole cone Cz. 



7 Discussion 

In this paper we have described relation between the classical MShG equation and its linear problem, 
on one hand, and quantum sine- and sinh-Gordon models on the other. This relation generalizes 
the relation |201I26| between ordinary differential equations [241125] and integrable structures of 
Conformal Field Theories j21fl23j to the massive case. We believe it also brings useful insight into 
the emergence of TBA equations in recent analysis of the MShG equation jl4fll8j . 

The discussion in this paper is in terms of the Q- and T-functions of the quantum models (II. 2p 
and ()1.13p . which are defined through the Bethe Ansatz for the vacuum states. More generally, these 
functions have to be understood as the vacuum eigenvalues one-parameter families of commuting 
operators Q(^) and T{6). In integrable lattice models (e.g. XXZ and XYZ chains) these operators 
were discovered in pioneering works of Baxter |32p33j. When an integrable quantum field theory 
emerges as continuous (scaling) limit of such lattice systems, it inherits these operators. However, 
for many reasons (including subtleties of the continuous limit) it is desirable to have constructions 
of these operators directly in field theoretic terms. There was some progress in this direction. Thus, 
in Ref. [21] the operators Tj(6'), j = |, 1, |, . . . were constructed explicitly for massless (conformal) 
field theory, as traces of quantum monodromy matrices over 2j + 1 dimensional auxiliary spaces. 
This construction of Tj{9) admits more or less direct extension to the massive sine-Gordon model 
and its reductions [41J- In both massless and massive cases it allows one to establish directly basic 
properties of these operators, which we summarize in Appendix B. In particular, one can argue 
that all Tj(0) are entire functions of 0, in the sense that all their simultaneous eigenvalues are entire 
functions. Furthermore, in the massless case similar construction (with auxiliary space supporting 
representation of g-oscillator algebra) exists for the Q-operator |221l23j. It is plausible that it also 
admits generalization to the massive case, but the details were never elaborated. However, both 
the massless construction and the lattice theory suggest certain properties of the sine-Gordon Q- 
operator. Thus, one expects that Q(^) is entire function of 9 as well. This operator commutes with 
all Tj{9'), and satisfies the famous T — Q equation of Baxter, 

Ti(0)Q(0) = Q(0+f)+Q(e + f) , (7.1) 

familiar from the lattice theory |32y33j . This equation is finite difference analog of a second order 
differential equation. Since in the sine-Gordon model T{9) is periodic function of 9 (see Eq. ()B.7p ). 
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one expects to have two "Bloch wave" solutions Q±(0), 



Q+(0 + M^)=UQ+(0), (7.2) 

with some unitary operator U commuting with Q±{9). Again, by comparison to the massless 
limit [22l[23llll], it is natural to identify U with the Flouquet-Bloch operator associated with the 
discrete symmetry of the sine-Gordon theory, 

Vipix,t)V-^ = ipix,t) + 2Tr/^ . (7.3) 

Since Ti(^) is invariant with respect to the charge conjugation (IB.lh . the operators Q+(0) and 

2 

Q_(0) are related to each other by this symmetry transformation, Eq. ()B.13P (this symmetry was 
already taken into account in writing (17.2p ). so one can deal with one independent operator, say 
Q = Q-I-. Additional piece of analytic information - the asymptotic behavior 

log Q{9) ~ const e^^ + 0(1) (7.4) 

as ?R.e{9) — )• ±oo in the strips H± (I3.10j) - can be inferred from the massless Q-operators j22|[23]. 
These assumptions and some of their simple consequences (in particular, precise form of (|7.4p ) are 
summarized in Appendix B. 

Once the sine-Gordon Q-function Q{9,k) is understood as the fc-vacuum eigenvalue of the 
Q-operator, the question arises about its eigenvalues associated with the excited states. Basic 
properties of such eigenvalues can be inferred from ()B.10p - (]B.15p . Can the excited-state eigenvalues 
be also related to integrable classical equations? In the massless case it is known how to modify the 
Schrodinger equation (I2.3ip to accommodate for the excited states [31j. In the massive theory this 
is interesting open question. Of course, generalization of the DDV equation ()3.32p to the excited 
states is well known [TOlllSj. 
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Appendix 

A Derivation of Eqs.(l3l])-(|33D 

Here we sketch derivation of the properties (I3.ip - (I3.4|) of Q±. 
Eqs.(l32|) fohow from ([2T9|) . (f2:2U]) and easily established relation 

a^S*{p,(j)\-9) = -S{p,(l)\9) at Sm(0) = . (A.l) 
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To prove the quasiperiodicity (j3.ip we use the evident relation 



Q-{e) = W det a{p,(t)\e), 
where W = — cos(7r/). Then 



det 



i7r(Q+l) ' 

20L , 



W X 



(A.2) 



(A.3) 



The following formula is immediate consequence of Eqs. (|2.16|) . (|2.17|) . 

(p, <A + i I e - ) = e-' (p, <A I 9) 

It is also straightforward to prove similar relation for the solution H, 

Using (IXill and (|X5l) one has 



-ie''"' 1^ X 



(A.4) 



(A.5) 



(A.6) 



det 



2a 



le 



i7r(a+l) 
2a 

Q-{e + 



)• 



2a 



+ 



i7r(o+l) 
2a 



i(ci+l)7r 
2a 



For given values of p and the solution 3 (|2.23p . considered as the function of complex 6, 
is analytic in the strip — ^ — (a + 1) </> < Qm{9) < ^ — (a + 1)0. Since is entire function 
of 9, one concludes that Q+{9) is analytic in the strip 13=771(0)1 < vr + and hence due to the 
quasiperiodicity (|3.ip it is entire function of 9. Since Q± and 'if± are entire functions of 9, H 
()2.25p is an entire function as well. 

To prove the Quantum Wronskian relation (j3.3p we use (j2.25p to obtain 



Q- 



2a' 



+ Q+i9±i^)^.{p,^Ti^\9±^). 

From (j2.16p we have 

S{p, T ^ I ± f ) = [ Q„(0 ± f ) <A I 0) + Q+{e ± ^) cp 1 1 

and hence 

det[s{p,cl>-^\e + ^),S{p,cP + ^\0-^)] = det 

{Q-{9 + 1) Q+iG -t)-Q-i(^- 1) Q+i(^ + t)) ■ 

Using normalization condition (12.18P and easily established relation 



(A.7) 



(A.8) 



det 



one arrives at (j3.3p . 

Finally, Eqs. ()3.4p are immediate consequences of (j2.2ip . 



-2i , 



(A.9) 



(A.IO) 



31 



B T and Q -operators in sine-Gordon model 

B.l T-operators 

Let C and P be unitary operators of charge conjugation and parity transformation in the sine- 
Gordon model ()1.2p . 

Cip{x,t)C = -(p{x,t) , (B.l) 
Fip{x,t)F = (p{-x,t) . (B.2) 

Integrabihty of the quantum sine-Gordon model can be expressed in terms of family of operators 
("transfer-matrices") Tj{d) , 6* S C , j = ^, 1, |, • • •, having the following properties (a = — 1): 



Mutual commutativity 



U and C invariance 



Parity transformation 



Hermiticity 



Periodicity 



[Tj{e),T,,{0')] = O . (B.3) 



[Tj{e),U] = [T,{e),C] = . (B.4) 



'Tj{e)F = Tj{-e) . (B.5) 



T](0) = Tj(r). (B.6) 



Tj(0+i^^) =Tj(0) . (B.7) 

Fusion relation 

T.(0)T,(^+M|±1))=T^._,(0+M|±H))+T^.^,(0 + |H) . (B.8) 
Asymptotic at real 9: 

logTi(0)~^ 2(-l) sm(^2^)C„x<j _ ^ « _ " (^'9) 



n=0 



2n-i ^ as y — )• — oo 



Here I_i = I_i = ^ is a c-number, while l2n-i (n = 1, 2. . .) are local IM (j3.59p . and the 
constants C„ are given by (|3.63p . 

The operators Tj{0) can be thought of as the continuous limits of the Baxter's commuting transfer- 
matrices [33j, or traces of monodromy matrices of the quantum inverse scattering method [34j . 
Construction of these operators directly in continuous quantum field theory is outlined in Refs. 

[HHiaiii]. 
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B.2 Q-operators 



Here we summarize expected properties of the operators Q±{0). Eqs. ()B.10p . (|B.lip just recapitulate 
what was already suggested in Section[71 The rest follows from the T—Q equation (j7.ip . the defining 
relations (|7.2p . and the asymptotic ()7.4p . 



Commutativity 



(B.IO) 



mvariance 



Quantum Wronskian relation 

where U is the field translation (I7.3p . 
Charge and parity conjugations 

CQ±(0)C = Q^(0) , 



2a J 



u 



Hermiticity 



Leading asymptotic in the strips ff+ (j3.1ip (we assume a 7^ 1, 3, 5 . . .) 



U=^2 §2 exp 



U=^t S^t exp 



. 4cos(^) _ 
Mi2e-^±^ 



as 5Re (0) — > +00 , 
as (6*) — )■ —00 



(B.ll) 



(B.12) 



(B.13) 



(B.14) 



(B.15) 



4cos(^) 

Here S is some operator which is invariant with respect to the P and U symmetries, 

[S, P] = [S, U] = , (B.16) 

and satisfies the relations 

S-^=CSC, S"f = S. (B.17) 

At the moment we do not know physical interpretation of the operator S, but regard it as 
very interesting open question. In the massless case it is similar to the Liouville "reflection S- 
matrix" [47^1 . 



We note in this connection that in the massless case the fuU spectrum of S can be extracted from recent remarkable 
paper [48] . 
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